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Abstract 

According to very recent developments of the LARES mission, which would be devoted to the mea- 
surement of the general relativistic Lense-Thirring effect in the gravitational field of the Earth with 
Satellite Laser Ranging, it seems that the LARES satellite might be finally launched in a polar, low- 
altitude orbit by means of a relatively low-cost rocket. The observable would be the node only. The 
Lense-Thirring effect on it would consist of a secular linear trend. The biasing classical secular nodal 
precessions due to the even zonal harmonics of the geopotential, which represent the major source 
of uncertainty, exactly vanished if and only if the orbit would be exactly polar. Due to the small 
altitude, even small possible deviations from the projected inclination, which might be induced by the 
orbital injection errors, would yield a rather large systematic error due to the even zonal harmonics of 
geopotential in the measurement of the relativistic nodal shift. So, in this paper we show how such a 
configuration, in fact, to the present level of knowledge of the terrestrial gravitational field according 
to the EGM96 gravity model, should be of relatively little utility in increasing the obtainable accuracy 
in measuring the Lense-Thirring effect with respect not only to the originally proposed supplementary 
LARES-LAGEOS configuration, but also to the present LAGEOS-LAGEOS II experiment which has 
a total accuracy of the order of 20% — 30%. Maybe the situation might be improved, at least to a 
certain extent, when the new, more accurate Earth gravity models from the CHAMP and GRACE 
missions will be available. 
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1 Introduction 

In its weak-field and slow-motion approximation General Relativity predicts that, among other 
things, the orbit of a test particle freely falling in the gravitational field of a central rotating 
body is affected by the so called gravitomagnetic dragging of the inertial frames or Lense- 
Thirring effect. More precisely, the longitude of the ascending node Q and the argument of the 
perigee u of the orbit undergo tiny precessions according to [Lense and Thirring, 1918] 

^LT = ^^'^ 3 , (1) 

c^a^i^l — e^) 2 

GGJcosi . , 

'^LT = f) (2) 

c^a^(l — e^) 2 

in which G is the Newtonian gravitational constant, J is the proper angular momentum of 
the central body, c is the speed of light in vacuum, a, e and i are the semimajor axis, the 
eccentricity and the inclination, respectively, of the orbit of the test particle. The Lense- 
Thirring precessions for the LAGEOS satellites amount to 

fiLAGEOs ^ 3imas/y, (3) 

^LAGEOS II ^ 3i.5inas/y, (4) 

^LAGEOS _ 31.6mas/y, (5) 

^LAGEOS II ^ _57nias/y. (6) 

The first measurement of this effect in the gravitational field of the Earth has been obtained 
by analyzing a suitable combination of the laser-ranged data to the existing passive geodetic 
satellites LAGEOS and LAGEOS II [Ciufolini et al, 1998]. The observable [CiufoUm, 1996] 
is a linear trend with a slope of 60.2 milliarcseconds per year (mas/y in the following) and 
includes the residuals of the nodes of LAGEOS and LAGEOS II and the perigee of LAGEOS 
II[]. The total relative accuracy of the measurement of the solve-for parameter /ilt; introduced 
in order to account for this general relativistic effect, is of the order of 2 x 10^^ [Ciufolini et 
al, 1998]. 



^The perigee of LAGEOS was not used because it introduces large observational errors due to the smallness 
of the LAGEOS eccentricity [Ciufolini, 1996] which amounts to 0.0045. 
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In this kind of satellite-based space experiments the major source of systematic errors is 
represented by the aliasing trends due to the classical secular precessions [Kaula, 1966] of the 
node and the perigee induced by the mismodelled even zonal harmonics of the geopotential 
J2, J4, Jq, ■■■ Indeed, according to the present knowledge of the Earth's gravity field based on 
the EGM96 model [Lemoine et al., 1998], they amount to a large part of the gravitomagnetic 
precessions of interest, especially for the first two even zonal harmonics. In the performed 
LAGEOS-LAGEOS II Lense-Thirring experiment the adopted observable allows for the can- 
cellation of the static and dynamical effects of J2 and J4. The remaining higher degree even 
zonal harmonics affects the measurement at a 13% level. 

In order to achieve a few percent accuracy, in [CiufoUni, 1986] it was proposed to launch 
a passive geodetic laser-ranged satellite- the former LAGEOS III - with the same orbital pa- 
rameters of LAGEOS apart from its inclination which should be supplementary to that of 
LAGEOS. 

This orbital configuration would be able to cancel out exactly the classical nodal precessions, 
which are proportional to cosi, provided that the observable to be adopted is the sum of the 
residuals of the nodal precessions of LAGEOS III and LAGEOS 

6n^^^ + 60} = 62/XLT. (7) 

Later on the concept of the mission shghtly changed. The area-to-mass ratio of LAGEOS III 
was reduced in order to make less relevant the impact of the non-gravitational perturbations 

and the eccentricity was enhanced in order to be able to perform other general rclativistic tests: 
the LARES was born [CiufoUni, 1998]. The orbital parameters of LAGEOS, LAGEOS II and 
LARES are in Table 1. 

Table 1: Orbital parameters of LAGEOS. LAGEOS 11 LARES and POL ARES. 



Orbital parameter 


LAGEOS 


LAGEOS 11 


LARES 


POLARES 


a (km) 


12,270 


12,163 


12,270 


8,378 


e 


0.0045 


0.014 


0.04 


0.04 


i (deg) 


110 


52.65 


70 


90 



Recent developments of the concept of the twin satellites in supplementary orbits have led 
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to the discover of new, possible gravitomagnetic observables based on the use of the perigees as 
well [lorio, 2002a] and of unexpected connections with the gravitomagnetic clock effect [lorio 
and Lichtenegger, 2002]. 

Unfortunately, at present we do not know if the LARES mission will be approved by any 
space agency. Although much cheaper than other proposed and/or approved space-based mis- 
sions, funding is the major obstacle in implementing the LARES project. The most expensive 
part is the launching segment. 

Very recently the possibility of launching the LARES satellite into an orbit with a = 8, 378 
km, i = 90 deg, e = 0.04 and using as observable its node has been considered [Lucchesi 
and Paolozzi, 2001]. In the following we will name POLARES the LARES satellite in such 
proposed polar orbit. The Lense-Thirring secular rate of the POLARES node would amount 
to 96.9 mas/y. Then, the observable would be 

Sif^ = 96.9/xlt, (8) 

i.e. a linear secular trend with a slope of 96.9 mas/y. The choice of a so low altitude is motivated 
by the need of using a cheap rocket launcher; so, it must certainly be considered as an admirable 
and important further effort towards the practical realization of the LARES project. The polar 
orbit would allow to prevent the ahasing effects of the mismodelled classical secular precessions 
of the node induced by the even / = 2n zonal m = coefficients of the multipolar expansion of 
the static part of the geopotential. The non-gravitational perturbations [Lucchesi, 2001; 2002; 
Lucchesi and Paolozzi, 2001] would represent a minor problem. 

In this paper we wish to critically analyze such important and interesting evolution of the 
LARES concept. 

2 The POLARES 

2.1 The node— only scenario 

Here we will show that the proposed POLARES configuration and the use of the node only as 
gravitomagnetic observable would present some drawbacks due to the gravitational static and 
time-varying perturbations. 
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In regard to the classical secular precessions induced by the centrifugal oblateness of the 
Earth, it should be considered that the choice of the inclination would be of crucial importance. 
Such orbital element is affected neither by the secular perturbations induced by the even zonal 
harmonics of the geopotential [Kaula, 1966] nor by the semisecular 18.6-year and 9.3-year tidal 
perturbations [lorio, 2001]. On the other hand, due to possible orbital injection errors which 
are closely related to the quality of the rocket launcher to be used, the POLARES inclination 
would be far from being exactly 90 deg. The relatively low altitude of POLARES would enhance 
the impact of even small departures of zpL from its projected nominal value on the systematic 
error due to the classical static even zonal harmonics of geopotential 6fiff^^ zonais_ fjj-^jg jg dearly 
shown in Figure 1 in which it has been calculated by means of the covariance matrix of the 
EGM96 Earth gravity model up to degree / = 20. It should be considered that, contrary to the 

5n'''-=96.9n : sensitivity to orbital injection errors in inclination (EGM96 up to 1=20) 





constant combination coefficient=1 ; ap|^=8,378 km; ep|^=0.04 






POLARES inclination (deg) 

Fi gure 1; POLARES scenario: influence of the orbital injection errors in the POLARES inclination on the 
zonal error in the Lense-Thirring nodal shift. 

LAGEOS satellites which are almost insensitive to the even zonal harmonics of degree higher 
than / = 20, this would not be the case for the POLARES, due to its projected altitude of 
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2,000 km. Then, it might turn out that the estimates of Figure 1 are optimistic. 

If we consider that the impact of the mismodelled even zonal harmonics of the geopotential 
on the current LAGEOS-LAGEOS II Lense-Thirring experiment is of the order of 13%, ac- 
cording to EGM96, a certain weakness of the POLARES node-only scenario becomes apparent. 
In regard to the possibilities offered by the forthcoming new gravity models from the CHAMP 
and GRACE missions, it should be considered that the major improvements in the accuracy 
of the even zonal coefficients of the geopotential are expected for the first degrees. Now, in 
the LAGEOS-LAGEOS II Lense-Thirring experiment the first two even zonal harmonics are 
cancelled out and the terms of degree higher than / = 20 are not relevant because the LAGEOS 
and LAGEOS II orbits are not sensitive to them. So, the systematic error due to the geopoten- 
tial, in this case, is due to the even zonal harmonics from / = 6 to / = 20. Then, the CHAMP 
and GRACE missions would certainly improve such error. On the other hand, the POLARES 
configuration would be sensitive both to the first two even zonal harmonics, which should be 
greatly improved by the CHAMP and GRACE results, and to the terms of degree higher than 
/ = 20, for which the improvements should be less relevant. So, even in this perspective, maybe 
that a simple reanalysis of the LAGEOS and LAGEOS II data should be more fruitful than 
the node-only POLARES choice. 

Such conclusion is enforced also by a simple evaluation of the impact of some tidal pertur- 
bations0 [/ono, 2001]. Putting aside the fact that the 18. 6- year lunar tide, which has a period 
of 18.6 years due to the lunar node only motion and a large amplitude, is a / = 2, m = 
constituent which would affect the POLARES node for ipL 7^ 90 deg, let us draw our attention 
to the 1 = 2, m = 1 Ki tesseral tide. It is an important constituent which exerts a relevant 
perturbing action on the node of a satellite^ and has the same period of just the node of the 
satellite. It turns out that for, say, ipL = 89.8 deg the period of the POLARES node would 
amount to 26,801 days, i.e. 73.3 years. This means that, even by not considering at all the even 
zonal perturbations, the mismodelled part of the Ki tide itself would resemble a superimposed 
aliasing trend over a reasonable observational time span of a few years and would completely 

^The amplitude of the tidal perturbations on the node is proportional to {1/ sin i)dFimp/di where Fimp{i) 
are the inclination functions [Kaula, 1966]. For the even zonal (m = 0), tesseral (m = 1) and sectorial (m = 2) 
terms they are F201 = (3/4)sin^i — (1/2), F211 = —(3/2) sin z cos i and F221 = (3/2)sin^i, respectively. 

•^For LAGEOS it has a nominal amplitude of the order of 10'' mas [lorio, 2001] 
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mask the Lense-Thirring trend. 

2.2 The combined residuals scenario 

Let us try to see what could happen by inserting the node and the perigee of POLARES in a 
suitable combination of orbital residuals with the nodes of LAGEOS and LAGEOS II and the 
perigee of LAGEOS II along the lines sketched in [Ciufolini, 1996; lorio, 2002b]. Recall that if 
N orbital elements are present in such combinations, the effects of the first N — 1 even zonal 
harmonics of the geopotential are cancelled out, irrespectively of the orbital geometry of the 
employed satellites. 

It turns out that by combining the nodes of LAGEOS, LAGEOS II and POLARES and the 
perigees of LAGEOS II and POLARES in 

SCl^ + ci50" + C26Cf^ + C3(5cj" + C45uj^^ = XuriUjT (9) 

the inclination = 90 deg is singular in the sense that the slope of the relativistic effect 
diverges because the coefficient C2 with which the node of POLARES enters the combination 
diverges. For values of ipi, close to 90 deg Figure 2 and Figure 3 show that the systematic error 
due to the remaining 5Jio, SJu, according to the covariance matrix of EGM96 up to degree 
I = 20, would seriously affect the measurement of the Lense-Thirring linear trend. Moreover, 
the impact of all the gravitational and non-gravitational time-varying perturbations would be 
greatly enhanced by the large values of C2. The situation would not be better even if we would 
include only the perigee of POLARES in the combined residuals 

SCl^ + ci50" + csJcI;" + C36u^^ = Xlt/^lt, (10) 

as shown by Figure 4. 

3 Conclusions 

In this paper the proposal of putting the LARES satellite into a polar, elliptical orbit with an 
altitude of 2,000 km in order to look at its gravitomagnetic secular node shift has been critically 
analyzed. 
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Fi gure 2: POLARES, LAGEOS, LAGEOS II scenario: influence of the orbital injection errors in the PO- 
LARES inclination on the zonal error in the combined residuals with the node of POLARES. 

The key point is that the mismodelled classical nodal secular precessions induced by the 
even zonal coefficients of the multipolar expansion of the terrestrial gravitational field vanished 
if and only if the inclination of the satellite would be exactly 90 deg. Of course, mainly 
due to possible orbital injection errors, this could never happen. It turns out that the low 
altitude of the proposed orbital configuration, and the consequent high sensitivity to the higher 
even degree zonal terms of the geopotential, would greatly enhance the impact of even small 
departures of the real values of the inclination from the nominal value of 90 deg. For example, 
for just ipL = 90 ± 0.2 deg the systematic gravitational error would be of the order of 5%-10%, 
according to the EGM96 Earth gravity model up to / = 20, which is of the same order of 
magnitude of the present LAGEOS-LAGEOS II experiment (Its total error, including various 
systematic gravitational and non-gravitational perturbations, is of the order of 20%-30%). Of 
course, such a situation would be further made critical by the possible use of a low-cost launcher 
which, unavoidably, would induce not negligible orbital injection errors. Moreover, the tesseral 
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Sn + Cj5n + CjSn +0260) + c^Sto =X)i|_p sensitivity to orbital injection errors in inclination 




91 91.5 92 92.5 93 93.5 94 94.5 95 

POLARES inclination (deg) 

Slope of the gravitomagnetic trend: sensitivity to orbital injection errors in inclination 
4000 I 1 1 1 1 1 1 1 1 




91 91.5 92 92.5 93 93.5 94 94.5 95 

POLARES inclination (deg) 



Fi gure 3: POLARES, LAGEOS, LAGEOS II scenario: influence of the orbital injection errors in the PO- 
LARES inclination on the zonal error in the combined residuals with the node of POLARES. 

Ki tidal perturbation, which has the same period of the satelhte's node, would induce a secular 
aliasing trend over an observational time span of a few years because its period would amount 
to several tens of years for near polar orbits. 

Another important point is that the proposed POLARES would not yield substantial im- 
provements also in the context of the combined residuals approach which allows to cancel out 
the contribution of the first even zonal coefficients of the geopotential irrespectively of the 
inclination of the satellites. Indeed, it turns out that a combination including the nodes of LA- 
GEOS, LAGEOS II, POLARES and the perigees of LAGEOS II and POLARES would be not 
defined for = 90 deg because the coefficient weighing the node of the LARES would go to 
infinity. For very small deviations of zpl from such a critical value the systematic error induced 
by the remaining even zonal harmonics would amount to 10%-20%. Moreover, the coefficient 
with which the node of POLARES would enter the combination would be much more larger 
than unity and would greatly enhance the impact of the gravitational and non-gravitational 
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Sn + c^Sn +c^S(i> + CgSto =X|i|_-|-: sensitivity to orbital injection errors in inclination 

36 I 1 1 1 1 1 1 1 1 r 
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POLARES inclination (deg) 

Slope of ttie gravitomagnetic trend: sensitivity to orbital injection errors in inclination 
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53 95 ' ' ' ' ' ' ' ' ' ' ' 

89 89.2 89.4 89.6 89.8 90 90.2 90.4 90.6 90.8 91 

POLARES inclination (deg) 



Figure 4: POLARES, LAGEOS, LAGEOS II scenario: influence of the orbital injection errors in the PO- 
LARES inclination on the zonal error in the combined residuals without the node of POLARES. 

time-dependent perturbations. Even the use of the nodes of LAGEOS and LAGEOS II and the 
perigees of LAGEOS II and POLARES, in which case ipL = 90 deg would not create problems, 
would not yield benefits because the systematic gravitational error would be of the order of 
20%-30%. 

Moreover, the accuracy of the practical data reduction from such version of LARES would 
be affected by the atmospheric drag. 

Finally, we can conclude that the proposed low-cost version of the LARES mission would 
not yield any significant improvements in the measurement of the elusive Lense-Thirring effect, 
according to the present-day level of knowledge of the Earth's gravitational field. Perhaps, the 
situation could improve to a certain extent with the new, more accurate gravity models from 
CHAMP and GRACE missions which should become available in the next few years. On the 
contrary, the original concept of the couple of supplementary satellites would deserve grater 
attention thanks to its much richer spectrum of high accuracy relativistic observables. 
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